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It is shown that the most important source of energy in 
ordinary stars. is the reactions of carbon and nitrogen with 
protons. These reactions form a cycle in which the original 
nucleus is reproduced, wz. C?+H=N8, NB=C8-+ et, 
CB+H=N4 NYtH=08, OU=NE+4 et, NU4+H=C? 
-+He!. Thus carbon and nitrogen merely serve as catalysts 
for the combination of four protons (and two electrons) 
into an a-particle (§7). 

The carbon-nitrogen reactions are unique in their 
cyclical character (§8). For all nuclei lighter than carbon, 
reaction with protons will lead to the emission of an 
a-particle so that the original nucleus is permanently 
destroyed. For all nuclei heavier than fluorine, only 
radiative capture of the protons occurs, also destroying the 
original nucleus. Oxygen and fluorine reactions mostly lead 
back to nitrogen. Besides, these heavier nuclei react much 
more slowly than C and N and are therefore unimportant 
for the energy production. 

The agreement of the carbon-nitrogen reactions with 
observational data (§7, 9) is excellent. In order to give the 
correct energy evolution in the sun, the central tempera- 
ture of the sun would have to be 18.5 million degrees while 


§1. INTRODUCTION 


HE progress of nuclear physics in the last 
| few years makes it possible to decide rather 
definitely which processes can and which cannot 
occur in the interior of stars. Such decisions will 
be attempted in the present paper, the discussion 
being restricted primarily to main sequence 


stars. The results will be at variance with some 


current hypotheses. 

The first main result is that, under present 
conditions, no elements heavier than helium can 
be built up to any appreciable extent. Therefore 
we must assume that the heavier elements were 
built up before the stars reached their present 
state of temperature and density. No attempt 
will be made at speculations about this previous 
state of stellar matter. 

The energy production of stars is then due 
entirely to the combination of four protons and 
two electrons into an a-particle. This simplifies 
the discussion of stellar evolution inasmuch as 


* Awarded an A. Cressy Morrison Prize in 1938, by the 
New York Academy of Sciences. 


integration of the Eddington equations gives 19. For the 
brilliant star Y Cygni the corresponding. figures are 30 
and 32. This good agreement holds for all bright stars of 
the main sequence, but, of course, not for giants. 

For fainter stars, with lower central temperatures, the 
reaction H+H=D-+e?t and the reactions following it, are 
believed to be mainly responsible for the energy produc- 
tion. (§10) | 

It is shown further (§5-6) that no elements heavier than 
Het can be built up in ordinary stars. This is due to the fact, 
mentioned above, that all elements up to boron are disin- 
tegrated by proton bombardment (a-emission!) rather than 
built up (by radiative capture). The instability of Be 
reduces the formation of heavier elements still further. 
The production of neutrons in stars is likewise negligible. 
The heavier elements found in stars must therefore 
have existed already when the star was formed. 

Finally, the suggested mechanism of energy production 
is used to draw conclusions about astrophysical problems, 
such as the mass-luminosity relation (§10), the stability 
against temperature changes (§11), and stellar evolution 


(§12). 


the amount of heavy matter, and therefore the 
opacity, does not change with time. 

The combination of four protons and two 
electrons can occur essentially only in two ways. 
The first mechanism starts with the combination 
of two protons to form a deuteron with positron 
emission, vtz. 


H+H=D-+et. (1) 


The deuteron is then transformed into He‘ by 
further capture of protons; these captures occur 
very rapidly compared with process (1). The 
second mechanism uses carbon and nitrogen as 
catalysts, according to the chain reaction 


C24 H=N8+ +4, N8 = C8 et 
CH PESEN S ys 
N44H=0"+ y, 
N3+H=C®+Het, 


Ol = N!5 + et (2) 


The catalyst C!? is reproduced in all cases except 
about one in 10,000, therefore the abundance of 
carbon and nitrogen remains practically un- 
changed (in comparison with the change of the 
number of protons). The two reactions (1) and 
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(2) are about equally probable at a temperature 
of 16-108 degrees which is close to the central 
temperature of the sun (19-10° degrees!). At 
lower temperatures (1) will predominate, at 
higher temperatures, (2). 

No reaction other than (1) or (2) will give an 
appreciable contribution to the energy produc- 
tion at temperatures around 20-10° degrees such 
as are found in the interior of ordinary stars. 
The lighter elements (Li, Be, B) would “burn” 
in a very short time and are not replaced as is 
carbon in the cycle (2), whereas the heavier 
elements (O, F, etc.) react too slowly. Helium, 
which is abundant, does not react with protons 
because the product, Li>, does not exist; in fact, 
the energy evolution in stars can be used as a 
strong additional argument against the existence 
of He® and Li® (§3). 

Reaction (2) is sufficient to explain the energy 
production in very luminous stars of the main 
sequence as Y Cygni (although there are diffi- 
culties because of the quick exhaustion of the 
energy supply in such stars which would occur 
on any theory, §9). Neither of the reactions (1) 
or (2) is capable of accounting for the energy 
production in giants; if nuclear reactions are at 
all responsible for the energy production in these 
stars it seems that the only ones which could 
give sufficient energy are 


H?+H =He? (3) 
Li®: 7-+H = He?: 4+ He’. 


It seems, however, doubtful whether the energy 
production in giants is due to nuclear reactions 
at all.’ 

We shall first calculate the energy production 
by nuclear reactions (§2, 4). Then we shall prove 
the impossibility of building up heavier elements 
under existing conditions (§5—-6). Next we shall 
discuss the reactions available for energy pro- 
duction (§5, 7) and the results will be compared 
with available material on stellar temperatures 
and densities (§8, 9). Finally, we shall discuss the 
astrophysical problems of the mass-luminosity 
relation (§10), the stability of stars against 
temperature changes (§11) and stellar evolution 
($12). 


1 B. Strömgren, Ergebn. d. Exakt. Naturwiss. 16, 465 
(1937). 
2 G. Gamow, private communication. 
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§2. FORMULA FOR ENERGY PRODUCTION 


The probability of a nuclear reaction in a gas 
with a Maxwellian velocity distribution was first 
calculated by Atkinson and Houtermans.’ Re- 
cently, an improved formula was derived by 
Gamow and Teller.4 The total number of 
processes per gram per second ist 


4 PX1X2 Tr g i 
35/2 MıMo h 





Here p is the density of the gas, xıxz the concen- 
trations (by weight) of the two reacting types of 
nuclei, mm, their masses, Zie and Zæ their 
charges, mM=mme/(mi+m:2) the reduced mass, 
R the combined radius, 


a=h?/me’Z,Z>2 (5) 


the “Bohr radius”? for the system, T/A the 
probability of the nuclear reaction, in sec.~, 
after penetration, and 


memeZ eZ 9” i 
=3(-——) (6) 
2h°kT 


If we measure p in g/cm’, T in volts and 7 
in units of 10° degrees, we have 


P=5.3- 10" pxix2eP (Zi, Zo) re go'sec.~', (7) 
+ =42.7(Z;22)4(A/T)}, (8) 


1 8R\? 
o=——_(—} e2(8R/a)? (9) 
A;A2(Z12Z2A )8 a 


where A,A¢ are the atomic weights of the reacting 
nuclei (A;=m;/my), A=m/my, my=mass of 
hydrogen. For the combined radius of nuclei 1 
and 2 we put 


R=1.6-10-(A,+A.)# cm. (10) 


Then we obtain for ¢ the values given in Table I. 
The values of ¢ for isotopes of the same element 
differ only very slightly. 

The values of T for reactions giving particles 
can be deduced from the observed cross sections 


s R. d'E. Atkinson and F. G. Houtermans, Zeits. f. 
Physik 54, 656 (1929). 
: T ae and E. Teller, Phys. Rev. 53, 608 (1938), 
tq. (3). 
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TABLE I. Values of ¢ for various nuclear reactions. 





R 

REACTION (10738 cm) Y REACTION R p 
H? +H 2.3 0.38 | Si? 4-H 5.0 29.3 
H?:+H 2.55 0.48 | C7-4+-H 5.4 75 
He!+H 2.75 0.81 | H?+H? 2.55 0.67 
Li7+H 3.2 0.91 | Be7+H? 3.3 1.18 
Be? +H 3.45 1.16 | Be?7+He? 3.45 7.9 
BU+H 3.65 1.52 | He*-+H? 2.9 0.57 
C2+H 3.75 2.00 | He++ He? 3.0; 1.09 
N¥#+H 3.95 2.718 | Het+ Het 3.2 1.29 
Ol+H 4.1 3.80 | Li?-+He4 3.55 4.9 
Fl+H 4.35 5.5 | Be?’+Het 3.5; 13.2 
Ne2+H 4.5; 7.7 | Be8-++Het 3.6; 16.2 
Me?6+ H 4.8 13.2 | C#+Het 4.0 230 


TABLE II. Cross sections and widths for some nuclear reac- 
tions giving particles. 


E o R 


REACTION REF. KV CM? CM r 

H2 +H2 =He? +n: z 100 1.7. 10726 2.6. 10783 3» 105 

L+H! =2He4 42 1.7 < 107830 3.2. 107: 4-104 

Lif+H! =Het+He asia ~ same as Li?7++H in natural Li 5+ 105 

target 

Lis-+H? = eat 7 212 19-107 3.2-10-8 4.106 

Li?+H?2 =2Het-+n 7 212 5.5* 10726 3.3+ 19718 107 
Lis + Het z = 

Be? +H! = Bee ene Ta 212 1.46107 3.5.1073  1.7.107 
Li? +Heż 

Be? +H? =, Be’ +-H3 7a 212 5+ 107% 3.6. 1078 6+ 105 
Bel0-+ H! 

BU+H! =3Het 7b 212 610-8 = 3.7-10713 2+106 


of such reactions with the use of the formula 
(cf. reference 4, Eq. (2)) 


TR? Aj +A. 32R 
ttt al (2) 
2E As a 


2re? 





22s] (11) 


hy 


where £ is the absolute energy of the incident 
particle (particle 1). Table II gives the experi- 
mental results for some of the better investigated 
reactions. In each case, experiments with low 
energy particles were chosen in order to make 
the conditions as similar as possible to those in 
stars where the greatest number of nuclear 
reactions is due to particles of about 20 kilovolts 
energy. The cross sections were in each case 
calculated from the thick target yield with the 
help of the range-energy relation of Herb, 


| 937), Ladenburg and M. H. Kanner, Phys. Rev. 52, 911 
1 
`H. D. Doolittle, Phys. Rev. 49, 779 (1936). 
1 J. H. Williams, W.G. Shepherd and R. O. Haxby, Phys. 
Rev. 52, 390 (1937). 
7a J, H. Willams, R. O. Haxby and W. G. Shepherd, 
Phys. Rev. 52, 1031 (1937). 
H. Williams, W. H. Wells, J. T. Tate and E. J. 
Hill Phys. Rev. 51, 434 (1937). 
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Bellamy, Parkinson and Hudson.® The widths 
obtained (last column of Table II) are mostly 
between 3-10° and 2-107 ev, with the exception 
of the reaction Li’-+-H = 2He! which is known to 
be ‘‘improbable.’’® 
The y-ray widths T, can be obtained from 
observed resonance capture of protons. Table III 
gives the experimental results. Two of the older 
data were taken from Table XXXIX of reference 
10; all the others are from more recent experi- 
ments on proton" 4 and neutron!ë capture. 
Although the results of different investigators 
differ considerably (e.g., for Li7+H!=Be3, I is 
between 4 and 40 volts the latter value being 
more likely) they seem to lie generally between 
about 4 and 40 volts. Ordinarily, the width is 
somewhat larger for the more energetic y-rays, 
as is expected theoretically. A not too bad 
approximation to the experiments is obtained by 
using the theoretical formula for dipole radiation 
(reference 10, Eq. (711b)) with an oscillator 
strength of 1/50. This gives 
T,~0.1E,2, (12) 


TABLE III. y-ray widths of nuclear levels. 


-RAY ENERGY 


REACTION REFERENCE WIDTH (VOLTS) (MEV) 
Liit+Hi=Beit+y | { n a 17 
Bu+tH!=CP++¥ 11 0.6 12, 16 
C8-+Hi=NB++y 12 0.6 2 
CB-+Hi= N+ y 13, 14 30 4,8 
FMY4H! = Ne%+ 7 10 0.6, 8, 18 6 
C+ yl=CB+ + 15 <2.5 5 
Ol6-+- ny) = OM + y 15 <2.5 4 


8s D. B. Parkinson, R. G. Herb, J. C. Bellamy and C. M. 
Hudson, Phys. Rev. 52, 75 (1937). 

9 M. Goldhaber, Proc. Camb. Phil. Soc. 30, 560 (1934). 

10 Hf, A. Bethe, Rev. Mod. Phys. 9, 71 (1937). 

iW. A. Fowler, E. R. Gaerttner and C. C. Lauritsen, 
Phys. Rev. 53, 628 (1938). 

2 R, B. Roberts and N. P. Heydenburg, Phys. Rev. 53, 
374 (1938). 

183 P, I. Dee, S. C. Curran and V. Petrzilka, Nature 141, 
642 (1938). The y-rays from C give about equally as 
many counts as those from C! The efficiency of the counter 
for C383 y-rays is about twice that for C?2, because the cross 
section for production of Compton and pair electrons is 
smaller by a factor 2/3 while the range of these electrons 
is about 3 times longer. With an abundance of C! of about 
1 percent, the y-width for this nucleus becomes 50 times 
that of C®, I am indebted to Dr. Rose for these calculations. 

14M. E. Rose, Phys. Rev. 53, 844 (1938). 

16 O, R. Frisch, H. v. Halban and J. Koch, Nature 140, 
895 (1937); Danish Acad. Sci. 15, 10 (1937). 
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where E, is the y-ray energy in mMU (milli- 
mass-units), and I, the y-ray width in volts. 
For quadrupole radiation, theory gives about 


(12a) 


Formulae (12), and (12a) will be used in the 
calculations where experimental data are not 
available; they may, in any individual case, be 
in error by a factor 10 or more but such a factor 
is not of great importance compared with other 
uncertainties. 

It should be noted that quite generally radia- 
tive processes are rare compared with particle 
emission. According to the figures given in 
Tables II and III, the ratio of probabilities is 
10*-10° in favor of particle reactions. 

In a number of cases, the reaction of a nucleus 
A with a heavy particle (proton, alpha-) must 
compete with natural -radioactivity of A or 
with electron capture. In those cases where the 
lifetime of radioactive nuclei is not known 
experimentally, we use the Fermi theory. 
According to this theory, the decay constant for 
G-emission is!® | 


T,~5-10-E,/ (quadrupole). 


B=0.9-10-4f(W) |G|? sec.—. (13) 


The matrix element G is about unity for strongly 
allowed transitions, and 


1 3 9 
SW) =(W— n(=wm->m-— 
30 20 15 


+2W log {W+(W?—1)*}, (13a) 
where W is the maximum energy of the £- 
particle, including its rest mass, in units of me? 
(m=electron mass). 

The probability of electron capture is 


Bo=0.9-10-*x?N(h/mc)?W?|G|? sec.-!, (14) 
where W is the energy of the emitted neutrino 
in units of mc? and N the number of electrons 
per unit volume. If the hydrogen concen- 
tration is xm, we have (reference 1, p. 482) 
N=6-103p-4(1+xu) (p the density), and 

Bo=1.5-10-4 p(1+xH) W?|G|? sec.-1. | (14a) 


16 C, L. Critchfield and H. A. Bethe, Phys. Rev. 54, 248, 
862 (L) (1938). 
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§3. STABILITY OF UNKNOWN ISOTOPES 


For the discussion of nuclear reactions it is 
essential to know whether or not certain isotopes 
exist (such as Li‘, Li5, Be®, Be’, B, B®, C!°, etc.). 
The criterion for the existence of a nucleus is its 
energetic stability against spontaneous disinte- 
gration into heavy particles (emission of a 
neutron, proton or alpha-particle). Whenever a 
light nucleus is energetically unstable against 
heavy particle emission, its life will be a very 
small fraction of a second (usually ~10-*° sec.) 
even if the instability is slight (e.g., Be® will 
have a life of 10-7! sec. if it is by 50 kv heavier 


- than two a-particles!’). 


For the question of the lifetime of radioactive 
nuclei, it is also necessary to know the mass 
difference between isobars. Similar information 
is required for estimating the y-ray width in 
capture reactions (cf. Eqs. (12), (12a)). 


H? and He? 


The most recent determination! of the re- 
action energy in the reaction H?-+-H?=He?+7! 
yielded 3.29 Mev as compared with 3.98 in 
H?-+-H?=H?+H?!, With a mass difference of 
0.75 Mev between neutron and hydrogen atom,’ 
this makes He’ more stable than H? by 0.06 Mev. 
This would be in agreement with the experi- 
mental fact that no H? is present in natural 
hydrogen to more than 1 part in 10". Even if 
He’ should turn out to be heavier than H’, the 
difference cannot be greater than about 0.05 
Mev=0.1 mc? which would make the life of He’ 
exceedingly long (~2000 years at the center of 
the sun, (Eq. (14a)), 3000 years in the complete 
atom, on earth). | 


Ht and Lit 


As was first pointed out by Bothe and 
Gentner,? it is definitely possible that Ht is 
stable. Lit is, of course, less stable because of. 
the Coulomb repulsion between its three protons. 
If it is stable, Lit would be formed when He? 
captures a proton and would thus play an 
important role in stars (cf. §5). The only possible 
estimate of the stability seems to be a comparison 


17H. A. Bethe, Rev. Mod. Phys. 9, 167 (1937). 

18 T. W. Bonner, Phys. Rev. 52, 685 (1938). 

19 H, A. Bethe, Phys. Rev. 53, 313 (1938). 

20 W., Bothe and W. Gentner, Naturwiss. 24, 17 (1936). 
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of Li‘Li® to the analogous pair B’B®. Reasonable 
estimates? make B® just unstable (0.3-0.7 
mMU), while B® comes out at the limit of 
stability (binding energy between —0.4 and 
+0.4 mMU), i.e., B8 about 0.3 to 0.7 mMU 
more stable than B®. Li® (see below) is unstable 
by 1.4-1.8 mMU; if one assumes the same 
difference, Li‘ is found to be unstable by ~1 
mMU. However, this argument is very uncertain 
and the possibility of a stable Lit cannot be 
excluded at present. H* would, from a similar 
argument, turn out stable by 0.6 mMU. 


He’ and Li’ 

The instability of He® is shown directly by the 
experiments of Williams, Shepherd and Haxby” 
_on the reaction 


Li’+ H?= He®+ Het. (15) 


From a measurement of the range of the a- 
particles, the mass of He? is found (reference 23, 
Table 73, p. 373) to be 5.0137 whereas the 
combined mass of an a-particle plus a neutron 
would be only 4.003 86+ 1.008 93 = 5.012 8. Thus 
He? is unstable by 0.9 mMU (milli-mass-units) 
which is far outside the experimental error 
(about 0.1-0.2 mMU). It might be argued that 
the a-particle group observed in reaction (15) 
might correspond to an excited state of He’. 
However, this is extremely improbable because 
a nucleus of such a simple structure as He’ 
(a-particle plus neutron) should not have any 
low-lying excited levels. (This holds both in 
the a-particle and the Hartree model of nuclear 
structure.) Moreover, it would be difficult to 
explain why the a-particle group corresponding 
to the ground state of He® should not have 
been observed. 

The conclusion that the mass of 5.0137 found 
by W. S. and H. really corresponds to the ground 
state of He® is supported by considerations of 
mass defects. In fact, the instability of He’ was 

21H, A. Bethe, Phys. Rev. 54, 436, 955 (1938). 

27. H. Williams, W. G. Shepherd and R. O. Haxby, 
Phys. Rev. 51, 888 (1937). 

23 M.S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 
247 (1937). 

24 With the possible exception of a doublet structure of 
the ground state, similar to Li’. However, the doublet 
separation should probably be much smaller than in Li’? 
because of the loose binding of He®, and presumably both 
components of the doublet are already contained in the 


rather broad a-particle group observed by Williams, 
Shepherd and Haxby. 
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first predicted by Atkinson?® on the basis of such 
considerations. Considering analogous nuclei, 
consisting of a-particles plus one neutron, we 
find that the last neutron is bound with a binding 
energy of 5.3 mMU in C® and with only 1.8 
mMU in Be’. A binding energy of —0.9 mMU 
in He? fits very well into this series while a 
positive binding energy would not. 

If He® is unstable, this is a fortiori true of Li’ 
since the binding energies of these two nuclei 
should differ by just the Coulomb repulsion 
between proton and alpha in Li>. This repulsion 
will be about 0.6-1 mMU, so that Li® is unstable 
by 1.4-1.8 mMU. 

Thus all the nuclear evidence?** points to the 
nonexistence of both He® and Li5. Even if there 
were no such evidence, astrophysical data them- 


selves would force us to this conclusion, because 


at a temperature of 2-10’ degrees (central 
temperature of sun) the energy production from 
the combination of He* and H forming Li® would 
be of the order of 10'° ergs/g sec. (cf. §4), as 
against an observed energy production of 2 
ergs/g sec. Only the nonexistence of Li> prevents 
this enormous production of energy. 


Be’ 


This nucleus is certainly unstable, as can be 
shown by comparing it with the known nucleus 
He® from which it differs by the interchange of 
protons and neutrons. The Coulomb energy 
which is the only difference between the binding 
energies of the two nuclei can be calculated 
rather accurately.” The instability against dis- 
integration into Het+2H is between 1 and 
2.6 mM U. 


Be’ 


This nucleus has been observed by Roberts, 


Heydenburg and Locher.?8 It decays with a 


2 R.d’E. Atkinson, Phys. Rev. 48, 382 (1935). 

2a Note added in proof:—Recently, F. Joliot and I. 
Zlotowski (J. de phys. et rad. 9, 403 (1938)) reported the 
formation of stable He® from the reaction He*+H?=Heé 
+H!, The evidence is based upon the emission of singly 
charged particles of long range when heavy paraffin is 
bombarded by a-particles. However, the number of such 
particles observed was exceedingly small (only 6 out of a 
total of 126 tracks). Furthermore, the mass given for He® 
by Joliot and Zlotowski (5.0106) is irreconcilable with the 
stability (against neutron emission) of the well-known 
nucleus Hes, 

26 R, B. Roberts, N. P. Heydenburg and G. L. Locher, 
Phys. Rev. 53, 1016 (1938). 
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TABLE IV. Corrected and additional nuclear masses, 
and binding energtes. 


BINDING ENERGY 
NUCLEUS Mass (MMU) REFERENCE 
ni 1.008 93 19 
He: 3.016 99 5.87 18 
H! 4.025 4 0.6 +1 
He! 4.003 86 29 
Lit 4.026 9 —-1 =l 
Hed 5.013 7 —0.9 +0.2 23 
Lie 5.013 6 —1.6 +0.3 
Be® 6.021 9 —1.8 +0.8 21 
Be? 7.019 28 5.7 26 
Be® 8.007 80 — 0.08 +0.04 28 
B8 8.027 4 0.0 +0.4 21 
B? 9.0164 —0.5 +0.2 21 
C10 10.020 2 3.8 21 
N2 12.022 5 -24 3 0.0 0.9 21 
NW 13.010 08 2.03 19 
5.1 21 


ou 14.013 1 


half-life of 43 days (mean life 60 days) and 
probably only captures K electrons. Calculations 
of the Coulomb energy,” on the other hand, 
would make positron emission just possible 
(positron energy ~0.1 mMU). As a compromise, 
we assume that the mass of Be’ is just equal to 
Li’ plus two electrons, 1.e., 7.019 28. 


Be’ 


The instability of Be® against disintegration 
into two a-particles has been definitely estab- 
lished by Paneth and Glückauf?” who have shown 
that the Be® formed in the photoelectric disinte- 
gration of Be? disintegrates into 2 He‘. Kirchner 
and Neuert have confirmed this conclusion by 
investigating the products of the disintegration 
B!-+H = Be$!+Het. They found that frequently 
two a-particles enter the detecting apparatus 
simultaneously, with a small angle (less than 50°) 
between their respective directions of flight; this 
is just what should be expected if the Be® formed 
breaks up into two a-particles on its way to the 
detector. From the average angle between the 
two alphas, the disintegration energy of Bes 
(difference Be®— 2He*) was estimated as between 
40 and 120 kev.?8 


27 F, A. Paneth and E. Glückauf, Nature 139, 712 (1937). 

28 F, Kirchner and H. Neuert, Naturwiss. 25, 48 (1937). 

28a Note added in proof:—These conclusions are com- 
patible with the new measurements of S. K. Allison, 
E. R. Graves, L. S. Skaggs and N. M. Smith, Jr. (Phys. 
Rev. 55, 107 (1939)) on the reaction energy of Be®+H 
= Be’+ R. 

22K. T. Bainbridge, Phys. Rev. 53, 922(A) (1938). 
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Bs 
The existence is doubtful; calculation?! by 
comparison with its isobar Li® gives a binding 


energy between —0.4 and +0.4 mMU. This 


nucleus is not very important for astrophysics. 
B? 
B® is almost certainly unstable, as can be 


shown by calculating the difference in binding 
energy (Coulomb energy) between it and its 


isobar Be’. The theoretical instability is between 


0.3 and 0.7 mMU, 0.3 being almost certainly a 
lower limit. However, in view of the smallness 
of this instability, we shall at least discuss what 
would happen if B® were stable ($6). It will turn 
out that this would make almost no difference 
at “ordinary” temperatures (2-107 degrees) and 
not much even at higher ones (108 degrees). 
For these calculations, we shall assume B® to be 
stable with 0.2 mMU which seems generous. 
C10 

C! is stable with 4 mMU against Be+2H. 
N?2 

N? is doubtful, mainly because the binding 
energy of its isobar, B”, is known only very 
inaccurately (between 2 and 3.3 mM U). Assum- 
ing 2 mMU, N® would probably be instable, 
with 3.3 stable. 

Table IV summarizes the binding energies of 
doubtful nuclei, and also gives some nuclear 


masses supplementary to and correcting those 
given in reference 23, p. 373. 


§4. REACTION RATES AT 2-107 DEGREES 


We are now prepared to actually calculate the 
rate of nuclear reactions under the conditions 
prevailing in stars. We choose a temperature of 
twenty million degrees, close to the temperature 
at the center of the sun. In order to have a 
figure independent of density and chemical 
composition, we calculate (cf. 7) 


P= (me2/x2)P/px1, (16) 


Ppx, gives the probability (per second) that a 
given nucleus of kind 2 undergoes a reaction 
with any nucleus of kind 1. If there are no other 
reactions destroying or producing nuclei of 
kind 2, 1/Ppx; will be the mean life of nuclei of 
kind 2 in the star. 
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TABLE V. Probability of nuclear reactions at 2-10" degrees.** 


REACTION Q (MMU) T (Ev) 
H+H =H?+ e+ 1.53 Ref. 16 
H?-+-H = He? 5.9 1E 
H:-+H = Het 21.3 I0 E 

He!-+H = Li4* (0.5) 0.02 D 
He!+-H = Lid* (0.2) 0.005 D 
Li6-+H = He!+ He 4.1 5-1058 X 
Li7+H =2 He! 18.6 4-104X 
Be?+H = B®? (0.5) 0.02 D 
Be? +H =Lit+ Het 2.4 106 X 
B9+H = C!o* 3.5 2D 
BuUtH=CH 9.2 10 D 
BU+H =3 Het 9.4 10°F 
CUu+tH=N” (0.4) 0.02 D 
C:2-+H =N! 2.0 0.6 X 
cC8-+H =N" 8.2 30 X 
N4¥+H =O 7.8 5D 
N¥4+H =C-+ Het 5.2 10° & 
O-H =F! 0.5 0.02 D 
F19-+H = Ol6+ Hed 8.8 10°F 
Ne#-++ H = Na?8 10.7 10 D 
Mg?6-+ H = Al?? 8.0 10 D 
S130- H = Pê! 7.0 10 D 
Cl37-+ H = A38 12.0 10 D 
H?+H?=He?+n 3.5 3-108 X 
Be7+H?=B** 18.5 10 D’ 
Be?’+H?=B°+n* 11.9 106 E 
Be’?+He= C!° 16.2 1D’ 
H2+ Het = Li® 1.7 4-10-30 
He3-+ Het = Be’ 1.6 0.02 D’ 
He!-+ Het = Be’* (0.05) 5-107? Q 
Li?7+He!= BU 9.1 1D’ 
Be’+ He! =C! 8.0 1p’ 
C+ Het = 0!6 7.8 1 QO’ 


P (sEc.7}) LIFE, FOR px1 =30 
12.5 8.5.1072 1.2-10! yr. 
13.8 1.3 -1072 2 sec. 
14.3 1.7-107! 0.2 sec. 
22.7 3-107? 1 day 
23.2 6-1078 6 days 
31.1 7-1073 5 sec. 
31.3 6-1074 1 min. 
38.1 6:1073 2000 yr. 
38.1 4.10 15 min. 
44.6 2.1078 5000 yr. 
44.6 107-1? 1000 yr. 
44.6 1.2-1077 3 days 
50.6 10717 108 yr. 
50.6 4.10716 2.5-106 yr. 
50.6 2.1071 5-104 yr. 
56.3 2.10717 5-107 yr. 
56.3 5-107 2000 yr. 
61.6 8.1072 1022 yr. 
66.9 4.10717 3-107 yr. 
71.7 5.1072 2.108 yr. 
81.3 10-26 101 yr. 
90.4 4.10730 3-102 yr. 
103.1 5-107% 2.10% yr. 
15.7 103 
45.9 2.1078 
50.7 2.10710 
80.5 3.10728 
27.5 3.10" 
47.3 3.10717 3-107 yr. 
50.0 _ 107% 
71.0 2.5.1072 
86 3.10730 3-102 yr. 
119 7:1078 





** The letters in the column giving the level width mean: X =experimental value; D =calculated for dipole radiation, from Eq. (12); D’ = 
dipole radiation with small specific charge, 1/4 to 1/20 of Eq. (12); @ =quadrupole radiation, Eq. (12a); and E =estimate. 


* These reactions are not believed to occur since their product or one of the reactants is unstable. 


cussion. 


Table V gives the results of the calculations, 
based on Eqs. (7) to (9). In the first column, 
the nuclear reactions are listed. All reactions 
which seemed of importance in the interior of 
stars were considered ; in addition, some reactions 
with heavier elements (O'* to Cl8”) were included 
in order to show the manner in which the 
reaction rate decreases. Moreover, seven reac- 
tions were listed in spite of the fact that their 
products or reactants are believed to be (§3) 
unstable (starred) or doubtful (question mark) ; 
these reactions are included in order to discuss 
the consequences if they did occur. 

The second column gives the energy evolution 
Q in the reaction, calculated from the masses 
(reference 23, Table LX XIII, and this paper, 
Table IV). In the third column, the width I 
determining the reaction rate (cf. §2) is tabu- 
lated. Wherever possible, this was taken from 
experiments (Tables II and III) or from the 
“empirical formulae” (12), (12a) for the radia- 


They are listed merely for the sake of dis- 


tion width. For the radiative combination of 
two nuclei of equal specific charge (H?+He+%, 
He*+ Het, C?-+He*) quadrupole radiation was 
assumed, otherwise dipole radiation.*° For almost 
equal specific charge (e.g., Be’-+Het*), the dipole 
formula with an appropriate reduction was used. 
In some instances, the width was estimated 
by analogy (eg., NY+H=C-+ Het) or from 
approximate theoretical calculations (H?+H 
= He’),3! The way in which T was obtained was 
indicated by a letter in each instance. 

The fourth column contains 7, as calculated 
from (8), the fifth P from (16). The wide varia- 
tion of P is evident, also the smallness of P for 
a-particle as compared with proton reactions. 

30 If the combined initial nuclei and the final nucleus 
have the same parity (as may be the case, e.g., for O'6++H 
= F1"), it is still possible to have a dipole transition if only 
the incident particle has orbital momentum one. This does 
not materially affect its penetrability if R>a (cf. (4), (5)) 


which is true in every case where the parities are expected 


to be the same. 
3 L, I. Schiff, Phys. Rev. 52, 242 (1937). 
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E.g., the reaction between Het and a nucleus as 
light as Be is as improbable as between a proton 
and Si. This arises, of course, from the greater 
charge and mass of the a-particle both of which 
factors reduce its penetrability. The reaction 
He!-+ He!= Be! has an exceedingly small proba- 
bility because of the small frequency and the 
quadrupole character of the emitted y-rays. 
Thus this reaction would not be important even 
if Be® were stable. On the other hand, the 
reaction Het+H = Li’ would be extremely prob- 
able if Li existed. The helium in the sun would 
be “burnt up” completely in about six days, 
even if rather unfavorable assumptions are made 
about the probability of the reaction. Similarly, 
if the energy evolution per process is 0.2 mMU 
= 3-10-7 ergs, the energy produced per gram of 
the star would be 


(6+ 10°8/4) px nxHe: 3+10-7-6- 10-8. 


With p=80, xH=0.35 and xye=0.1, this would 
give about 10" ergs/g sec. as against 2 ergs/g sec. 
observed. This is a very strong additional 
argument against the existence of Li?. 

In the last column of Table V, the mean life 
is calculated for the various nuclei reacting with 
protons, by assuming a density p=80 and 
hydrogen content xı=35 percent, which corre- 
spond to the values at the center of the sun.! 
It is seen that, with the exception of H, the 
lifetimes of all nuclei up to boron are quite short, 
ranging from a fraction of a second for H? to 
1000 years for B!°. (The life of B!° may actually 
be slightly shorter because of the reaction 
B'°+H = Be?+He*. See §6.) Of the two lives 
longer than 1000 years listed, one refers to B® 
which probably does not exist (§3), the other to 
Be’ which decays by positron emission with a 
half-life of 43 days.2 We must conclude that 
all the nuclei between H and C, notably H?, H3, 
Lif, Li’, Be’, BE, BY, can exist in the interior of 
stars only to the extent to which they are continu- 
ously re-formed by nuclear reactions. This con- 
clusion does not apply to He* because Li® does 
not exist. To He? it probably applies whether 
Lit exists or not, because He? will also be de- 
stroyed by combination with He‘ into Be’, 
although with a considerably longer period 
(3-10" years instead of the 1 day for the reaction 
giving Li‘). 
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The actual lifetime of carbon and nitrogen is 
much longer than it would appear from the table 
because these nuclei are reproduced by the 
nuclear reactions themselves (§7). This makes 
their actual lifetime of the order of 10” (or even 
102°, cf. §7) years, i.e., long compared with the 
age of the universe (~2-10° years). Protons, 
and all nuclei heavier than nitrogen, also have 
lives long compared with astronomical times. 


§5. THE REACTIONS FOLLOWING PROTON 
COMBINATION 


In the last section, it has been shown that all 
elements lighter than carbon, with the exception 
of H! and He!, have an exceedingly short life in 
the interior of stars. Such elements can therefore 
only be present to the extent to which they are 
continuously produced in nuclear reactions from 
elements of longer life. This is in accord with the 
small abundance of all these elements both in 
stars and on earth. 

Of the two more stable nuclei, He‘ is too inert 
to play an important rôle. It combines neither 
with a proton nor with another a-particle since 
the product would in both cases be an unstable 
nucleus. The only way in which Het can react at 
all, is by triple collisions. These will be discussed 
in the next section and will be shown to be very 
rare, as is to be expected. 

As the only primary reaction between ele- 
ments lighter than carbon, there remains there- 
fore the reaction between two protons, 


H!+H!=H?+ et. © 


According to Critchfield and Bethe," this process 
gives an energy evolution of 2.2 ergs/g sec. under 
“standard stellar conditions’ (2-107 degrees, 
p= 80, hydrogen content 35 percent). The reac- 
tion rate under these conditions is (cf. Table V) - 
2.5-10-19 sec.-!, corresponding to a mean life of 
1.2-10% years for the hydrogen in the sun. 
This lifetime is about 70 times the age of the 
universe as obtained from the red shift of 
nebulae. | 

According to the foregoing, any building up 
of elements out of hydrogen will have to start 
with reaction (1). The deuteron will capture 
another proton, 


H?+ H'= He’. (17) 
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This reaction follows almost instantaneously 
upon (1), with a delay of only 2 sec. (Table V). 
There is, therefore, always statistical equilibrium 
between protons and deuterons, the concentra- 
tions (by number of atoms) being in the ratio of 
the respective lifetimes. This makes the concen- 
tration of deuterons (by weight) equal to 


2-8.5.1072! 
x(H?) = x(H?!) ———— 
1.310 


=1.3-10-8x(H") (18) 


(cf. Table V). The relative probability of the 


reaction 





H?+H?=He?+ n! (19) 
as compared with (17), is then 
1 P(H?+H?=He'+n!) x(H?) 
4 P(H?4+H=He®) — x(H) 
108-1.3-107!8 
= 2.10, (19a) 
4-1.3-10-? 


(One factor 4 comes from the fact that in (19) 
two nuclei of the same kind interact; another is 
the atomic weight of H?.) Thus one neutron is 
produced for about 5-10" proton combinations. 
The further development of the He? produced 
according to (17) depends on the question of 
the stability of Lit and of the relative stability 
of H? and He. 


Assumption A: Li‘ stable 


In this case, the He’ will capture another proton, viz.: 
He3+ H!= Lif, (20) 


With the assumptions made in Table V, the mean life of 
He? would be 1 day. The Li‘ would then emit a positron: 


_ Lit=He!+er, (20a) 


With an assumed stability of Lit of 0.5 mMU compared 
with He®+H!, the maximum energy of the positrons in 
(20a) would be 20.8 mMU=19.4 Mev (including rest 
mass) which would be by far the highest 6-ray energy 
known. The lifetime of Lit may accordingly be expected 
to be a small fraction of a second (half-life = 1/500 sec. for 
an allowed transition in the Fermi theory). 

The most important consequence of the stability of Li‘ 
would be that only a fraction of the mass difference between 
four protons and an e-particle would appear as usable 
energy. For in the -emission (20a) the larger part of the 
energy is, on the average, given to the neutrino which will 
in general leave the star without giving up any of its energy 
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(see below). According to the Konopinski-Uhlenbeck 
theory, which is in good agreement with the observed 
energy distribution in 8-spectra, the neutrino receives on 
the average 5/8 of the total available energy if the latter 
is large. In our case, this would be 13.0 mMU. Adding 
0.2 mMU for the neutrino emitted in process (1), we find 
that altogether 13.2 mMU energy is lost to neutrinos, of a 
total of 28.7 mMU developed in the formation of an 
a-particle out of four protons and two electrons. Thus the 
observable energy evolution is only. 15.5 mMU, i.e., 54 
percent of the total. Therefore, if Li‘ is stable, process (1) 
would give only 1.2 ergs/g sec. instead of 2.2 (under 
“standard” conditions).16 

The neutrinos emitted will have some chance of pro- 
ducing neutrons in the outer layers of the star. It seems 
reasonable to assume that a neutrino has no other inter- 
action with matter than that implied in the 6-theory. Then 
a free neutrino (vy) will cause only “reverse 8-processes”?? 
of which the simplest and most probable is 


H-+v=n!+ et, (21) 


This process is endoergic with 1.9 mMU and is therefore 
caused only by fast neutrinos such as those from Li‘. The 
cross section is according to the Fermi theory 


o = 2"(h/mc)>-0.9-10-4c 1W(W2—1)3 
=1.7-10-8W(W2—1)* cm, (22) 


where W is the energy of the emitted positron, including 
rest mass, in units of mc®. In reaction (21), this is the 
neutrino energy minus 1.35 mMU. For the Li‘ neutrinos, 
the average cross section comes out to be 


ony = 2.5-107* cm? (22a) 


per proton, and the probability of process (21) for a neu- 
trino starting from the center of the star 


EE S R dr =1.5- 10-24 f P dr, (22b) 


where p(r) is the density (in g/cm?) at the distance 7 from 
the center of the star. For the sun, p=1.6-1077. This 
means that 1.6-1077 of the neutrinos emitted will cause 
reaction (21) before leaving the sun, and that the number 
of neutrons formed is 1.6-10~7 times the number of proton 
combinations (1). 

A further consequence of (20, 20a) would be that ordi- 
narily no nuclei heavier than 4 mass units are formed at 
all, even as intermediate products. Such nuclei would only 
be produced in the rare cases when H? or He? capture an 
a-particle rather than a proton, according to the reactions 

H? -+ Het = Lis (23) 
and 
He: + He! = Be’. (24) 
Under the favorable assumption that the concentration of 
He‘ is the same as of H! (by weight), the fraction of H? 
forming Li® is (cf. Table V) 


p(Li®) = 3-10710/1.3-107-7=2- 1078 (23a) 
the fraction of He® giving Be’ is 
p( Be’) =3-107 7/3 -10-7 = 10719, (24a) 


32 H. A. Bethe and R. Peierls, Nature 133, 689 (1934). 
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Most of the Li® will give rise to the well-known reaction 
Li‘+H = He’+ Het (23b) 


and most of the Be?’ will go over into Li’ which in turn 
reacts with a proton to give two a-particles. Only occasion- 
ally, Li®, Be? or Li? will capture an a-particle and thus 
form heavier. nuclei. It can be shown (cf. assumption B) 
that Li’ is the most efficient nucleus in this respect. 
Therefore, the amount of heavier elements formed is deter- 
mined by Be’, the mother substance of Li’, and is thus 107!° 
times the amount formed with assumption B. 


Assumption B: Lit unstable, He? more stable than H° 


This assumption seems to be the most likely according 
to available evidence. The only reaction which the He?’ 
can undergo, is then (24), i.e, each proton combination 
leads to the formation of a Be’? nucleus. The most probable 
mode of decay of this nucleus is by electron capture, leading 
to Li’. The lifetime of Be’ (half-life) is 43 days?® in the 
complete atom, and 10 months at the center of the sun 
(cf. 14, 14a). This makes the mean life=14 months and 
the reaction rate 2.8.1078 sec.. The capture of a proton 
by Be’ would, even if the product B8 is stable, be 2000 
times slower (Table V). Each electron capture by Be’ is 
accompanied by the emission of a neutrino of energy 
~2mc?=1.1 mMU (when Li’ is left in its excited state, 
which happens rather rarely, the neutrino receives only 
0.6 mMU). The total energy lost to neutrinos (including 
process 1) will therefore be very small in this case (~1.3 
mMU per a-particle formed, i.e., 4} percent of the total 
energy evolution) and practically the full mass energy will 
be transformed into heat radiation. The Li’ formed by 
electron capture of Be’ will cause the well-known reaction 


Li’+H =2 He’ (25) 


and have (Table V) a mean life of only 1 minute at 2-107 
degrees. 


The reaction chain described leads, as in the 
case of assumption A, to the building up of one 
a-particle out of four protons and two electrons, 
for each process (1). No nuclei heavier than 
He!‘ are formed permanently. Such nuclei can be 
produced only by branch reactions alternative to 
the main chain described. These will be discussed 
in the following. 


a. Reactions with protons.—When Li’ reacts with slow 
protons, the result is not always two a-particles, but, in 
one case out of about*® 5000, radiative capture, giving Be’. 
However, Be® will disintegrate again into two a-particles 
(§3), and during its life of about 107! sec., the probability 
of its reacting with another particle (e.g., capture of an- 
other proton) is exceedingly small (~10724}. Similarly, Be” 
will, in one out of about 2000 cases (see above) capture a 
proton and form B: if that ‘nucleus exists. However, B® 


33 It was assumed that radiative capture takes place only 
through the resonance level at 440 kv proton energy. The 
proton width of this level was taken as 11 kv, the radiation 
width as 40 ev. 
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will again go over into Be’, by positron emission, and two 
alphas will again be the final result. 

At this place, obviously, stability of Be® would increase 
the yield of heavy. nuclei. Then one stable Be® would be 
formed for 5000 proton combinations; and, if B® is also 
assumed to be stable, every Be® goes over into B®. Since 
about one out of 3-108 B® gives a C® (§6), the number of 
heavy nuclei (C!) formed would be ~107 per a-particle. 
This would be the highest yield obtainable. However, Be® 
is known to be unstable (§3). 

b. Reactions with a-particles—The only abundant light 
nucleus other than the proton is He*, The only reaction 
possible between an a-particle and Li’ or Be? is radiative 
capture, vz. 

Li’+Het=B", (26) 
Be?7+Het=C#, (26a) 


The probability of formation of B} and C" is (Table V) 


P(LIFHet) 2.54107. 
Biesen al, 26b 
p(B") P(L' +H) 6-10~ | oe 
p(CH) = P(Be?+He‘) 14 months = 4-101, (26c) 


~ P(Be?7+e=Li7) 3-10% yr. 


Thus the formation of B” is about as probable as that of 


CH; the effect of the lower potential barrier of Li’ for 


a-particles is compensated by its shorter life. The C” will, 
of course, also give B” by positron emission. 
The B! will react with protons in two ways, viz. 


B'-+-H!=3 Het, (27) 
BU+H!=C, (27a) 
The branching ratio is about 104: 1 in favor of (27) (calcu- 


lated from experimental data). Thus there will be one C!? 
nucleus formed for about 10% a-particles, The building up 


„of heavier nuclei, even in this most favorable case, is 


therefore exceedingly improbable. 

c. Reactions with He?,—Since He? has a rather long life 
3-107 years, Table V) and penetrates more easily through 
the potential barrier than the heavier He’, it may be con- 
sidered as an alternative possibility. However, the prob- 
ability of formation of C’ from Be’+ He is only 100 times 
greater than that of C! from Be’+Het* (Table V) if the 
concentrations of He® and He? are equal. Actually, that 
of He’ is only about 3-10~4 (life of He? divided by life of 
protons) so that this process is 1/30 as probable as (26a). 
For Li7+He3, the situation is even less favorable. 

d, Reactions with H?.—Deuteron capture by Be’ would 
lead to B® whose existence is very doubtful. The probability 
per second would be (cf. Table V and Eq. (18)) 


2-1078px(H?) =2-10~-30-1.3-10718=10779, (28) 


which is only 1/10 of the probability of (26a) (Table V). 
Moreover, most of the B® formed reverts to He‘ (§6) so 
that the contribution of this process is negligible. 

e. Reactions with Het in statu nascendt.—The process 
(25) produces continuously fast a-particles which need not 
penetrate through potential barriers. These a-particles 
have a range of 8 cm each in standard air, corresponding 
to 16 cm=0.02 g/cm? for both. In stars, with their large 
hydrogen content, a somewhat smaller figure must be used 
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since hydrogen has a greater stopping power per gram; 
we take 0.01 g/cm?. The cross section for fast particles is 
about 


(29) 


With r=1 volt (Table V) and R=3.6:10718 cm (Eq. (10)), 
this gives o=1.3-1073! cm?, The number of Be’ atoms per 
gram is 


14 months 
1.2-10" yr 


with x(H) =0.35, This gives for the number of processes 
(26a) per proton combination: 


0.01 -2.102.1.3-1073 = 2.5.1072, (29b) 


which is about the same as the formation of C” or B" by 
capture of slow alphas (cf. 26b, c). 

Returning to the main reaction chain in the case of our 
assumption B, we note that the formation of Be’ (Eq. 
(24)) is a very slow reaction, requiring 3-107 years at 
“standard” conditions (2-107 degrees). At lower tempera- 
tures, the reaction will be still slower and, finally, it will 
take longer than the past life of the universe (~2-10° 
years). In this case, the amount of He? present will be 
much smaller than its equilibrium value (provided there 
was no He? “in the beginning’) and the energy production 
due to reactions (24), (25) will be reduced accordingly. 
Ultimately, at very low temperatures (<12-10° degrees), 
the reaction H+H will lead only to He’, and will therefore 
give an energy production of only 7.2 mMU, i.e., only 
one-quarter of the high temperature value, 27.4 mMU. 


6-10%x(Be7) =6- 10%: x(H)=2-10% (29a) 


Assumption C: H: more stable than He? 
In this case, He® will be able to capture an electron, 
He®+ e = H8. 
Under the assumption that a difference in mass of 0.1 
electron mass exists between He* and Hi, the probability 

of (30) is, according to (14a), 

p(H3) =1.5-1071! sec.7! (30a) 
for a density p= 80 and 35 percent hydrogen content. This 
corresponds to a mean life of ~2000 years. The electron 


capture is therefore about 104 times more probable than 
the formation of Be’ according to (24). This ratio will be 


reversed at temperatures >4-107 since (30) is independent 


of T and the probability of (24) increases as T™. 
H°? will capture a proton and form Hef, 


H+H = Het 


with a mean life of about 0.2 seconds. This way of forma- 
tion of He! from reaction (1) is probably the most direct 
of all. As in B, practically no energy is lost to neutrinos. 

The formation of heavier elements goes as in B, but now 
there is only one Be’ formed for 104 proton combination 
processes. This reduces the probability of formation of C” 
by another factor 104, to one C! in 10?8 alphas. 

The H? itself does not contribute appreciably to the 
building up of elements. It is true that the reaction Be1"-+ H° 
=B%+7n! is about 100 times as probable as Be7+ H?= B® 


(30) ` 
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(considering the shorter life of H*), and therefore 10 times 
as probable as (26a). However, most of the B® reverts to 
He? (cf. §6) so that (26) and (26a) remain the most efficient 
processes for building up C¥. 


Summarizing, we find that the formation of 
nuclei heavier than Het can occur only in 
negligible amounts. One C! in 10% a-particles 
and one neutron in 10'* a-particles are the 
yields when Li‘ is unstable, one C! in 10* alphas 
and one neutron in 10’ alphas when Li‘ is stable. 
The reason for the small probability of formation 
of C® is twofold: First, any nonradioactive 
nucleus between He and C, ie., Li® 7, Be’, 
B® 1. reacting with protons will give a-particle 
emission rather than radiative capture so that a 
disintegration takes place rather than a building 
up. This will no longer be the case for heavier 
nuclei so that for these a building up is actually 
possible. Second, the instability of Be® causes a 
gap in the list of stable nuclei which is the harder 
to bridge because Be® is very easily formed in 
nuclear reactions (small mass excess). On the 
other hand, the instability of He® and Li® is of 
no influence because Be’ and Li’ are stages in 
the ordinary chain of nuclear reactions. 


§6. TRIPLE COLLISIONS OF ALPHA-PARTICLES 


In the preceding section, we have shown that 
collisions with protons alone lead practically 
always to the formation of a-particles. In order 
that heavier nuclei be formed, use must there- 
fore certainly be made of the a-particles them- 
selves. However, collisions of an a-particle with 
one other particle, proton or alpha, do not lead 
to stable nuclei. Therefore we must assume 
triple collisions, of which three types are con- 
ceivable: 


He!+2H = Be’, (31) 
2Het+H =B’, (32) 
3Het= C”, (33) 


The first of these reactions leads to a nucleus 
which is certainly unstable (Be). Even if it 
were stable, it would not offer any advantages 
over Be’ which is formed as a consequence of 
the proton combination (1). The second reaction 
leads to B® which is probably also unstable. 
However, since this is not absolutely certain, we 


34(5, Breit and E. Wigner, Phys. Rev. 49, 519 (1936). 
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shall discuss this process in the following. The 
last process leads directly to C', but since it 
involves a rather large potential barrier for the 
last a-particle, it is very improbable at 2-10’ 
degrees (see below). 


The formation of B® 


The probability of this process is enhanced by the well- 
known resonance level of Be’, which corresponds to a 
kinetic energy of about E=50 kev of two a-particles. The 
formation of B® occurs in two stages, 


2 He‘=Be’, Be’+H=B?2 (32a) 


with a time interval of about 107 sec, (life of Be’). The 
process can be treated with the usual formalism for 
resonance disintegrations, the compound nucleus being 
Bes, This nucleus can “disintegrate” in two ways, (a) into 
two a-particles, (b) with proton capture. We denote the 
respective widths of the Be’ level by 'g and I'y; the latter 
is given by the ordinary theory of thermonuclear processes, 
i.e., 


Ty =hpmoe/xe, (34) 


where ? is given by (4), and the subscripts 1 and 2 denote 
H?! and Be’, respectively. 
The cross section of the resonance disintegration becomes 
then 
Spaen Ol s 
(E-E) +} Ta tTn) 


"i (35) 
E, is the resonance energy. I'g is much larger than rH 
(corresponding to about 10" and 1071! sec,—, respectively) 
but very small compared with E, (about 107? against 5-104 
volts). The resonance is thus very sharp, and, integrating 
over the energy, we obtain for the total number of processes 


per cm’ per sec. simply 
pp=B(E,)v,7wA,2aT Hy. (36) 


Here B(E)dE is the number of pairs of a-particles with 
relative kinetic energy between E and E+dE per cm, viz. 


—E(kT 





Xa\? 2 Et 
) (36a) 


B(E) =3( Sant 


Me 


(X2=concentration of Het by weight). Combining (36), 
(36a) and (34), (4), we find 


1673 Pxaxa Trad 
352 malmy (kT)}R 
Xexp {4(2R/a)}— r— E,/kT}, 








p(B?) = 


aR?r? 
(37) 


where Traq is the radiation width for the process Be84-H 
=B°*, Numerically, (37) gives for the decay constant of 
hydrogen the value 


pmy/xH = 1.00-1074(pxe)*P eT fe 7B,/T eT, (31a) 


where T is measured in millions of degrees, the resonance 
energy E, of Be® in kilo-electron-volts and T in ev. The 
quantities T, and 7 refer to the process Be8+H = B°, If 
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there were no resonance level of Be’, (37) would be replaced 
by 
8 p*x.?xH Trad 
— ——— —— aR? ' R”? 2/2 
243 mma h ONTT 
Xexp {(32R/a)?+ (32R'/a')3—r—r'}, (38) 
where the primed quantities refer to the reaction 2 He‘ 
= Be’, the unprimed ones to Be8+H=B°. Numerically, 
(38) gives 
pmya/xu = 2.0107 (px%a)T ore To'r eT., (38a) 


Assuming T=20, p=80, xa=0.25, T=0.02 electron- 
volts, we obtain for the probability of formation of B? per 
proton per second: 


p(B’) = 


pmn/xu=2.1073 for resonance, Er=25 kev 
10-*! for resonance, E, = 50 kev 
4.1078 for resonance, E,=75 kev 
2:107% for resonance, Z,=100 kev 
5-10- for nonresonance. 


The value 25 kev for the resonance level must probably be 
excluded on the basis of the experiments of Kirchner and 
Neuert.?8 But even for this low value of the resonance 
energy, the probability of formation of B® is only 107° 
times that of the proton combination H+H=H?+ et 
(Table V, px: =30). With Z,=50 kev which seems a likely 
value, the ratio becomes 4-107!5, On the other hand, the 
building up of B® (if this nucleus exists) would still be the 
most efficient process for obtaining heavier elements (see 
below). 


Reactions of B® 


It can easily be seen”! that B® cannot be positron-active 
but can only capture electrons if it exists at all. If B® is 
stable by 0.3 mMU, the energy evolution in electron 
capture would be just one electron mass. The decay con- 
stant of B® (for 8-capture) is then, according to (14a), 
1.5-10-9 sec. (9 =80, xy =0.35) corresponding to a life- 
time of about 20 years. On the other hand, the lifetime 
with respect to proton capture (Table V) is 5000 years. 
Therefore, ordinarily B® will go over into Be®. This nucleus, 
in turn, will in general undergo one of the two well-known 
reactions: 

Be®+H = Be’+-H?, (39) 
Be®-++-H = Li§+ He‘. (39a) 


Only in one out of about 10° cases, B!° will be formed by 
radiative proton capture. Therefore the more efficient way 
for building up heavier elements will be the direct proton 
capture by B°, leading to C!°, which occurs in one out of 
about 300 cases, 

The C!® produced will go over into B!® by positron 
emission. BY may react in either of the following ways: 


B09+H = 01, (40) 
Bl0-++ H = Be’+ Het, (40a) 


The reaction energy of (40a) is (cf. Table VII, §8) 1.2 
mMU; the penetrability of the outgoing alphas about 1/40 
(same table), therefore the probability of the particle reac- 
tion (40a) will be about 100 times that of the capture 
reaction (40). 
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The CH from (40) will emit another positron. The result- 
ing B} reacts with protons as follows: 


B!+H =C”, (41) 
BU+H =3 Hel, (41a) 


Both reactions are well-known experimentally. Reaction 
(41) has a resonance at 160 kev. From the width of this 
resonance and the experimental yields, the probability of 
(41) with low energy protons is about 1 in 10,000 (i.e., the 
same as for nonresonance). Altogether, about one B® in 
3-108 will transform into C®, 

With a resonance energy of Be’ of 50 kev, and 2-10’ 
degrees, there will thus be about one C” formed for 10?! 
a-particles if B® is stable. This is better than any other 
process but still negligibly small. 

At higher temperatures, the formation of B® will become 
more probable and will, for T> 108, exceed the probability 
of the proton combination. At these temperatures (actually 
already for T>3-107) the B°?’ will rather capture a proton 
(giving C!°), Even then, there remain the unfavorable 
branching ratios in reactions (40), (40a) and (41), (41a),® 
so that there will still be only one C!2 formed in 10° alphas. 
Thus even with B® stable and granting the excessively high 
temperature, the amount of heavy nuclei formed is ex- 
tremely small. 


Direct formation of C” 


C2 may be formed directly in a collision between 3 
a-particles. The calculation of the probability is exactly 
the same as tor the formation of B®. The nonresonance 
process gives about the same probability as a resonance 
of Be® at 50 kev. With p =80, x,=i, '=0.1 electron-volt, 
T =2-107 degrees, the probability is 1075 per e-particle, 
i.e., about 10787 of the proton combination reaction (1). 
This gives an even smaller yield of C than the chains 
described in this and the preceding section. The process is 
strongly temperature-dependent, but it requires tempera- 
tures of ~10° degrees to make it as probable as the proton 
combination (1). 


The considerations of the last two sections 
show that there is no way in which nuclei heavier 
than helium can be produced permanently in 
the interior of stars under present conditions. 
We can therefore drop the discussion of the 
building up of elements entirely and can confine 
ourselves to the energy production which is, in 
fact, the only observable process in stars. 


§7. THE CARBON-NITROGEN GROUP 


In contrast to lighter nuclei, C is not perma- 
nently destroyed when it reacts with protons; 


% The reaction C'+H=N¥® becomes more probable 
than C= BU-+ et only at T>3-108 degrees. The branching 
ratio in (40), (40a) may perhaps be slightly more favorable 
because the effect of the potential barrier in (40a) may be 
stronger. 
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instead the following chain of reactions occurs 


(Ce4HI=NB, (42) 
N8 =CB4¢t, (42a) 
CB+HI=N", (42b) 
NU+Ht=0%8, (42c) 
OS = NB 4 ¢t, (42d) 
N4+H!=C2+ Het, (42e) 


Thus the C” nucleus is reproduced. The reason 
is that the alternative reactions producing a-par- 
ticles, vz. 


C2+H!=B?+Het, (43) 
C8-++H!= B+ He’, (43a) 
N¥4+H!=C"-+ Het, (43b) 


are all strongly forbidden energetically (Table 
VII, §8). This in turn is due to the much greater 
stability of the nuclei in the carbon-nitrogen 
group as compared with the beryllium-boron 
group, and is in contrast to the reactions of Li, 
Be and B with protons which all lead to emission 
of a-particles. 

The cyclical nature of the chain (42) means 
that practically no carbon will be consumed. 
Only in about 1 out of 104 cases, N‘® will capture 
a proton rather than react according to (42e). In 
this case, O'* is formed: 


N44+H!=01, (44) 


However, even then the C” is not permanently 
destroyed, because except in about one out of 
5-107 cases, O18 will again return to C® (cf. §8). 
Thus there is less than one C! permanently con- 
sumed for 10" protons. Since the concentration 
of carbon and nitrogen, according to the evidence 
from stellar spectra, is certainly greater than 
107! this concentration does not change notice- 
ably during the evolution of a star. Carbon and 
niirogen are true catalysts; what really takes place 
is the combination of four protons and two electrons 
into an a-pariicle. 

A given C! nucleus will, at the center of the 
sun, capture a proton once in 2.5:10° years 
(Table V), a given N** once in 5-107 years. 
These times are short compared with the age of 
the sun. Therefore the cycle (42) will have re- 
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TABLE VI. Central temperatures necessary for giving ob- 
served energy production in sun, with various nuclear 


reactions. 
T 
REACTION (MILLION DEGREES) 
H2+H = He? 0.36 
He!+H =Lié 2.1 
Li? +H =2He! 2.2 
Be?+H = Lit+Het 3.3 
BYe-+H =CH 9.2 
BU+H =3He! 5.5 
C24+H = N8 15.5 
N4+H =0% 18.3 
Ot -+H =F" 32 
Ne”-+H = Na?8 37 


peated itself many times in the history of the 
sun, so that statistical equilibrium has been 
established between all the nuclei occurring in 
the cycle, viz. CPCBNBN¥4N1O15.. In statistical 
equilibrium, the concentration of each nucleus 
is proportional to its lifetime. Therefore N" 
should have the greatest concentration, C! less, 
and C8N!5 still less. (The concentration of the 
radioactive nuclei N} and O!’ is, of course, very 
small, about 107! of N'*). A comparison of the 
observed abundances of C and N is not very 
conclusive, because of the very different chemical 
properties. However, a comparison of the isotopes 
of each element should be significant. 

In this respect, the result for the carbon 
isotopes is quite satisfactory. C captures slow 
protons about 70 times as easily as C” (experi- 
mental value!), therefore C” should be 70 times 
as abundant. The actual abundance ratio is 
94:1. The same fact can be expressed in a 
more “experimental” way: In equilibrium, the 
number of reactions (42) per second should be 
the same as of (42b). Therefore, if a natural 
sample of terrestrial carbon (which is presumed 
to reproduce the solar equilibrium) is bom- 
barded with protons, equally many captures 
should occur due to each carbon isotope. This is 
what is actually found experimentally; the 
equality of the y-ray intensities from C! and 
C8 is, therefore, not accidental.*$ 

The greater abundance of C! is thus due to 
the smaller probability of proton capture which 
in turn appears to be due to the smaller hy of 

36 It would be tempting to ascribe similar significance 
to the equality of intensity of the two a-groups from 
natural Li bombarded by protons (Li®+-H =He!+ He’, 


Li7-++-H =2He‘). However, the lithium isotopes do not 
seem to be genetically related as are those of carbon. 
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the capture y-ray. Thus the great energetic 
stability of C! actually makes this nucleus 
abundant. However, it is not because of a Boltz- 
mann factor as has been believed in the past, but 
rather because of the small energy evolution of 
the proton capture reaction. | 

In nitrogen, the situation is different. Here 
N" is énergetically less stable (has higher mass 
excess) than N?® but is more abundant in spite 
of it (abundance ratio ~500 : 1). This must be 
due to the fact that N** can give a p—a reaction 
while N” can only capture a proton; particle 
reactions are always much more probable than 
radiative capture. Thus the greater abundance of 
N" is due not to its own small mass excess but 
to the large mass excess of C” which would be 
the product of the p—a reaction (43b). 

Quantitative data on the nitrogen reactions 
(42c), (42e) are not available, the figures in 
Table V are merely estimates. If our theory 


about the abundance of the nitrogen isotopes is 


correct, the ratio of the reaction rates should be 
500 : 1, i.e., either N44+H=O must be more 
probable*** or N" +H =C!2+He* less probable 


than assumed in Table V. Experimental investi- 


gations would be desirable. 

Turning now to the energy evolution, we notice 
that the cycle (42) contains two radioactive 
processes (N and O") giving positrons of 1.3; 
and 1.8; mMU maximum energy, respectively. 
If we assume again that 2 of the energy is, on the 
average, given to neutrinos, this makes 2.0 mMU 
neutrino energy per cycle, which is 7 percent of 
the total energy evolution (28.7 mMU). There 
are therefore 4.0-10-5 ergs available from each 
cycle. (It may be mentioned that the neutrinos 
emitted have too low energy to cause the trans- 
formation of protons into neutrons according 
to (21).) 

The duration of one cycle (42) is equal to the 
sum of the lifetimes of all nuclei concerned, i.e., 
practically to the lifetime of N. Thus each N” 
nucleus will produce 4.0-10-5 erg every 5-107 
years, or 3-10-®° erg per second. Under the 


36a Note added in proof:—lIn this case, the life of N” in 
the sun might actually be shorter, and its abundance 
smaller than that of C£. Professor Russell pointed out to 
me that this would be in better agreement with the evidence 
from stellar spectra. Another consequence would be that a 
smaller abundance of N! would be needed to explain the 
observed energy production. 
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TABLE VII. p-a reactions. 


ENERGY POTENTIAL 
INITIAL PRODUCT EVOLUTION BARRIER PENETRABILITY 
NucLeus NUCLEUS Q (MMU) B (MMU) P 

Lié He? 4.14 1 

Li? Hef 18.59 1 
Be? Lié 2.45 2.6 1 

Bio Be? 1.2 3.4; 0.027 
Bu 2Het 9.4 E 
Cr B°? —8.1 0 

C13 Bio — 4.4 0 
N“ Cu — 3.5 0 
Nb CL 5.2 4.6 1 

Ole NB — 5.8 0 

QO Ni 1.3 5.1 1.7-1074 
Ol N15 3.1 5.0 0.2 
F19 Q16 8.8 1 
Ne? F1!” — 4.5 0 
Ne?! Fis — 1.6 0 
Ne? Fis —1.6 0 
Na? Ne?’ 1.6 6.7 6-1075 
Mg? Na” <— 3.0 0 
Mg” Na”? — 2.1 0 
Mg?6 Naz — 2.0 0 
Al?? Mg? 1.8 7.6 1.1-107° 
Si28 Al <—2.9 0 
Si? Ales — 2.0 0 
Si? AP? —2.4 0 
pst $128 2.0 8.5 2.5-1077 
S82 p29 <0 0 

S33 p30 sD 0 

S34 pP»: — 2.0 0 
Cl% s82 2.3 9.3 5-1078 
C137 5% 4.2 9.1 2.5-1078 


assumption of a N concentration of 10 percent 
by weight, this gives an energy evolution of 


6: 103.0.1 


a -3-10-29= 100 ergs/g sec. 


(45) 


at ‘‘standard stellar conditions,” i.e., T=2-10’, 
p= 80, hydrogen concentration 35 percent. 

This result is just about what is necessary to 
explain the observed luminosity of the sun. 
Since the nitrogen reaction depends strongly on 
the temperature (as 7'8) and the temperature, 
as well as the density, decrease rapidly from the 
center of the sun outwards, the average energy 
production will be only a fraction, perhaps 1/10 
to 1/20, of the production at the center.?®® This 
means that the average production is 5 to 10 
ergs/g sec., in excellent agreement with the ob- 
served luminosity of 2 ergs/g sec. 


366 Added in proof:—According to calculations of R. 
Marshak, the correct figure is about 1/30. 
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Thus we see that the reaction between nitrogen 
and protons which we have recognized as the 
logical reaction for energy production from the 
point of view of nuclear physics, also agrees 
perfectly with the observed energy production 
in the sun. This result can be viewed from 
another angle: We may ignore, for the moment, 
all our nuclear considerations and ask simply 


which nucleus will give us the right energy 


evolution in the sun? Or conversely: Given an 
energy evolution of 20 ergs/g sec. at the center 
of the sun, which nuclear reaction will give us 
the right central temperature (~19- 10° degrees) ? 

This calculation has been carried out in Table 
VI. It has been assumed that the density is 80, 
the hydrogen concentration 35 percent and the 
concentration of the other reactant 10 percent 
by weight. The “widths” were assumed the same 
as in Table V. Given are the necessary tempera- 
tures for an energy production of 20 ergs/g sec. 
It is seen that all nuclei up to boron require 
extremely low temperatures in order not to give 
too much energy production ; these temperatures 
(<10’ degrees) are quite irreconcilable with the 
equations of hydrostatic and radiation equi- 
librium. On the other hand, oxygen and neon 
would require much too high temperatures. 
Only carbon and nitrogen require nearly, and 
nitrogen in fact exactly, the central temperature 
obtained from the Eddington integrations (19-106 
degrees). Thus from stellar data alone we could 
have predicted that the capture of protons by 
N! is the process responsible for the energy 
production. : 


§8. REACTIONS WITH HEAVIER NUCLEI 


Mainly for the sake of completeness, we shall 
discuss briefly the reactions of nuclei heavier 
than nitrogen. For the energy production, these 
reactions are obviously of no importance because 
the higher potential barrier of the heavier nuclei 
makes their reactions much less probable than 
those of the carbon-nitrogen group. 

The most important point for a qualitative 
discussion is the question whether a p—a reac- 
tion is energetically possible for a particular nu- 
cleus, and, if possible, whether it is impeded by 
the potential barrier. In Table VII are listed the 
energy evolution in —a reactions for all stable 
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(nonradioactive) nuclei up to chlorine. In the 
first column, the reacting nucleus is given, in 
the second, the product of a p—a reaction. 
The third column contains the reaction energy Q; 
when Q is negative, the reaction is energetically 
impossible so that the initial nucleus can only 
capture protons with y-emission. In the fourth 
column, the height of the nuclear potential 
barrier is given for all reactions with positive Q. 
In the last column, the penetrability of the 
potential barrier is calculated according to 
standard methods (reference 10, p. 166). If 
Q>B, the penetrability is 1; if Q is negative, 
P=0 was inserted. 

The a priori probability of a —a reaction is 
roughly 10* times that of radiative capture. 
Therefore the emission of a-particles will be 
preferred when P>10-*. It is seen from the 
table that for all nuclei up to boron the p—a 
reaction is strongly preferred, a fact which we 
recognized as the main reason for the impossi- 
bility of building up heavier elements than 
He!’ (§5). Furthermore, in the carbon group, 
only proton capture is possible for C¥C8N'4 
while for N' the p—a reaction will strongly 
predominate (cf. §7). 

The oxygen-fluorine group shows intermediate 
behavior. O! can only capture protons, for O" 
the capture and the a-emission will have roughly 
equal probability while for O'8 and F!® the p—a 
reaction will be much more probable. With a 
ratio 10* for the a priert probabilities, about 40 
percent of the O" will become F'!8 (and then O18 
by positron emission) while 60 percent will 
revert to N'4. Of the O'8, only 1 part in 2000 will 
become F", and of the F!’, only 1 in 10,000 will 
transform into Ne?. Thus, under continued 
proton bombardment, about one O!* nucleus in 
5-10? will ultimately transform into Ne?®, the 
rest will become nitrogen. 

Actually, these considerations are somewhat 
academic because in general the supply of protons 
will be exhausted long before all the O!‘ initially 
present in the star will have captured a proton. 
Because of the small energy evolution in the 
reaction O'’4-H=F"’, this reaction is extremely 
slow (~10" years) so that equilibrium in the 
oxygen group will not be reached in astronomical 
times. 
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Among the nuclei heavier than fluorine, the 
p—a reaction is in general energetically per- 
mitted only for those with mass number 4n +3. 
But even for these, the energy evolution is so 
much less than the height of the barrier that the 
penetrabilities are extremely small. Thus for all 
these elements only proton capture will occur 
(with the possible exception of Cl”). 

These considerations demonstrate the unique- 
ness of the carbon-nitrogen cycle. 


§9. AGREEMENT WITH OBSERVATIONS 


In Table VIII, we have made a comparison 
of our theory (carbon-nitrogen reaction) with 
observational data for five stars for which such 
data are given by Strémgren.! The first five 
columns are taken from his table, the last con- 
tains the necessary central temperature to give 
the correct energy evolution with the carbon- 
nitrogen reactions (cf. Table VI). As in §7, we 
have assumed a N* content of 10 percent, and an 
energy production at the center of ten times the 
average energy production (listed in the second 
column). 

The result is highly satisfactory : The tempera- 
tures necessary to give the correct energy 
evolution (last column) agree very closely with 
the temperatures obtained from the Eddington 
integration (second last column). The only excep- 
tion from this agreement is the giant Capella: 
This is not surprising because this star has greater 
luminosity than the sun at smaller density and 
temperature; such a behavior cannot possibly be 
explained by the same mechanism which ac- 


TABLE VIII. Comparison of the carbon-nitrogen reaction with 
observations. 


CENTRAL TEM- 
PERATURE 
(MILLION DEGREES) 


CONTENT ENERGY 
LUMINOSITY CENTRAL  (PER- INTE- PRODUC- 
STAR ERG/GSEC. DENSITY CENT) | GRATION TION 
Sun 2.0 76 35 19 18.5 
Sirius A 30 4i 35 26 22 
Capella 50 0.16 35 6 32 
U Ophiuchi 180 12 50 25 26 
(bright) 
Y Cygni 1200 6.5 80 32 30 
(bright) 


37 A, S. Eddington, The Internal Constitution of the Stars 
(Cambridge University Press, 1926). 


450 H. A. 


counts for the main sequence. We shall come 
back to the problem of energy production in 
giants at the end of this section. 

For the main sequence we observe that the 
small increase of central temperature from the 
sun to Y Cygni (19 to 33. 10° degrees) is sufficient 
to explain the much greater energy production 
(104 times) of the latter. The reason for this is, 
of course, the strong temperature dependence of 
our reactions (~T8, cf. §10). We may say that 
astrophysical data themselves would demand such a 
strong dependence even if we did not know that 
the source of energy are nuclear reactions. The 
small deviations in Table VIII can, of course, 
easily be attributed to fluctuations in the nitro- 
gen content, opacity, etc. 

In judging the agreement obtained, it should 
be noted that the “observational” data in Table 
VIII were obtained by integration of an Edding- 
ton model,! 3” i.e., the energy production was 
assumed to be almost constant throughout the 
star. Since our processes are strongly tempera- 
ture dependent, the ‘“‘point source” model should 
be a much better approximation. However, it 
seems that the results of the two models are not 
very different so that the Eddington model may 
suffice until accurate integrations with the point 
source model are available.” 

Since our theory gives a definite mechanism of 
energy production, it permits decisions on ques- 
tions which have been left unanswered by 
astrophysicists for lack of such a mechanism. 
The first is the question of the ‘‘model,”’ which is 
answered in favor of one approximating a point 
source model. The second is the problem of 
chemical composition. The equilibrium condi- 
tions permit for the sun a hydrogen content of 
either 35 or 99.5 percent when there is no helium, 
and intermediate values when there is helium. 
The central temperature varies from 19-10° to 
9.5-10° when the hydrogen content increases 


87a Mr, Marshak has kindly calculated the central tem- 

perature and density of the sun for the point source model, 
using Strémgren’s tables for which we are indebted to 
Professor Strömgren. With an average atomic weight 
w=1, Marshak finds 

for the point source model 7, =20.3-105, p.=50.2 

for the Eddington model T,=19.6-108, p.=72.2 
Not only is the temperature difference very small (34 
percent) but it is, for the sake of the energy production, 
almost compensated by a density difference in the opposite 
direction. The product peT +8 is only 20 percent greater for 
the point source model. 
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from 35 to 99.5 percent. It is obvious that the 
latter value can be definitely excluded on the 
basis of our theory: The energy production due 
to the carbon-nitrogen reaction would be re- 
duced by a factor of about 108 (100 for nitrogen 
concentration, 10° for temperature). The proton 
combination (1) would still supply about 5 
percent of the observed luminosity; but apart 
from the fact that a factor 20 is missing, the 
proton combination does not depend sufficiently 
on temperature to explain the larger energy 
production in brighter stars of the main sequence. 
Thus it seems that only a small range of hydrogen 
concentrations around 35 percent is permitted; 
what this range is, depends to some extent on 
the N! concentration and also requires a more 
accurate determination of the distribution of 
temperature and density. 

Next, we want to point out a rather well-known 
difficulty about the energy production of very 
heavy stars such as Y Cygni. With an energy 
production of 1200 ergs/g sec., and an available 
energy of 1.0-10-5 erg per proton (formation of 
a-particles!), all the energy will be consumed in 
1.7-108 years. Since at present Y Cygni still has a 
hydrogen content of 80 percent, its past life 
should be less than 3.5-10’ years. We must 
therefore conclude that Y Cygni and other heavy 
stars were “born’’ comparatively recently—by 
what process, we cannot say. This difficulty, 
however, is not peculiar to our theory of stellar 
energy production but is inherent in the well- 
founded assumption that nuclear reactions are 
responsible for the energy production.*8 

Finally, we want to come back to the problem 
of stars outside the main sequence. The white 
dwarfs presumably offer no great difficulty. The 
internal temperature of these stars is probably 
rather low,! 3% because of the low opacity 
(degeneracy !) sothat the small energy production 
may be understandable. Quantitative calcu- 
lations are, of course, necessary. For the giants, 
on the other hand, it seems to be rather difficult 
to account for the large energy production by 
nuclear reactions. If the Eddington (or the point 


38 Even if the most stable nuclei (Fe, etc.) are formed 
rather than He, the possible life will only increase by 30 
percent. 
39S, Chandrasekhar, Monthly Not. 95, 207, 226, 676 

1935). 
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TABLE [X. Relations between stellar constants. 


-a_a 


y=18 y =18 y =15 y =4.5 
a=% a=1 a= a =l 
B=234 8 =3)4 B =234 B =234 
QUANTITY 6 = 19384 6 =2014 6 = 1634 6=6\4 
Radius R M? -Tut 5T ( yg)? .05 M? Tip? 51 (yg) 0 .05 M?.70,,0 49 (yz) 0 .06 M? 20-0 .36(+y2) 0.16 
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source) model is used, the central temperatures 
and densities are exceedingly low, e.g., for 
Capella T=6-10°, p=0.16. Only a nuclear re- 
action going at very low temperature is there- 
fore at all possible; Li7+-H=2He* would be 
just sufficient. But it seems difficult to conceive 
how the Li’ should have originated in all the 
giants in the first place, and why it was not 
burned up long ago. The only other source of 
energy known is gravitation, which would require 
a core model*® for giants.*4 However, any core 
model seems to give small rather than large 
stellar radii. 


§10. THE Mass-LuMINOSITY RELATION 


In this section, we shall use our theory of 
energy production to derive the relation between 
mass and luminosity of a star. For this purpose, 
we shall employ the well-known homology re- 
lations (reference 1, p. 492). This is justified 
because we assume that all stars have the same 
mechanism of energy evolution and therefore 
follow the same model. Further, it is assumed 
that the matter throughout the star is non- 
degenerate which seems to be true for all stars 
except the white dwarfs. (For all considerations 
in this and the following two sections, cf. 
reference 1.) 

We shall consider the mass of the star M, the 
mean molecular weight u, the concentration of 
“Russell Mixture” y and the product of the 
concentrations of hydrogen and nitrogen, z, as 
independent variables. In addition, we introduce 
for the moment the radius R which, however, 
will be eliminated later. Then, obviously, we have 
for the density (at each point) 


p~M/R. 


40, Landau, Nature 141, 333 (1938). 
41 This suggestion was made by Gamow in a letter to 
the author. 
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From the equation of hydrostatic equilibrium 
dp/dr= —GM,p/r* (47) 
(G= gravitational constant) and the gas equation 
b= RT p/u (47a) 


(R the gas constant, radiation pressure neglected), 
we find 


T~ Mu/R. (48) 


Finally, we must use the equation of radiative 
equilibrium : 


i Sk L; 
qgT'— = —-— 
dr 4c 4rr 





(49) 


where a is the Stefan-Boltzmann constant, c the 
velocity of light, k the opacity, and L, the 
luminosity (energy flux). at distance r from the 
center. 

For the opacity, we assume 


k~ yp T- (50) 


(y concentration of heavy elements). Usually, a 
is taken as 1 and 6=3.5 (Kramers’ formula). 
However, the Kramers formula must be divided 
by the “guillotine factor” + which was calculated 
from quantum mechanics by Strömgren. For 
densities between 10 and 100, and temperatures 
between 10’ and’ 3-10’, Strémgren’s numerical 
results can be fairly well represented by taking 


t~piT-. Therefore we adopt a=4, B=2.75 
in (50). | 
The luminosity may be written 


L~ MpzTY. (51) 


That the energy production per unit mass, L/M, 
contains a factor p follows from the fact that it is 
due to two-body nuclear reactions; this factor is 


42 Cf. reference 1, Table 6, p. 485. 
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Fic. 1. The energy production in ergs/g sec. due to the 
proton-proton combination (curve H+H) and the carbon- 
nitrogen cycle (N+H), as a function of the central tem- 
perature of the star. Solid curve: total energy production 
caused by both reactions. The following assumptions were 
made: central density =100, hydrogen concentration 35 
percent, nitrogen 10 percent; average energy production 
1/5 of central production for H+H, 1/10 for N+-H. 


apparent from all our formulae, e.g., (4). 
Z=% 1X2 is the product of the concentrations of the 
reacting nuclei (N“ and H). For y we obtain 
from (4), (6) 


d log (r’e~7) 
jG) 


(52) 
d log T 


For N“+H and T=2.-10', this gives y= 18. For 
T=3.2-107 (Y Cygni), y=15.5; generally, 
y~T-}. y=18 will be a good approximation over 
most of the main sequence. 


Inserting (50), (51) in (49), we have 
Tit- wyg M pte R. (53) 


Combining this with (46), (48), and introducing 
the abbreviation 


b= ytitsa- 8 (54) 
we find 
R ~M Ota) ld yl Tisa) S (yg) E, (55) 
T ~ M2 2+a) l? y O30) (yg) LS, (56) 
p ~ MH 2+a)/ô y3+3 (1+3) 18 (yg) -318 (57) 
L ~ M3tat2 (2+a) (8—30) /8 4+3 a+ (743a) (B—8a)/5 
xy Hye) HOI, (58) 


Furthermore, the surface temperature may be of 
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interest, we have 
T ş ~Li R~? ~ Mitat} (2+a)(8—8a+2)/8 


x< uzti att (743 a) (B—Ba+2)/day—% 


XK (yz) ee, (59) 


In Table IX, these formulas are given explicitly, 
for four different sets of constants. 

The most important result is that the central 
temperature depends only slightly on the mass of 
the star, viz. as M®™?5 and M®3° for y=18 and 15. 
The reason for this is the strong temperature 
dependence of the reaction rate : The exponent of 
M in (56) is inversely proportional to 6 which is 
mainly determined (cf. (54)) by the exponent y 
in formula (51) for the temperature dependence 
of the reaction rate. The integration of the 
Eddington equations with the use of observed 
luminosities, radii, etc., gives, in fact, only a 
small dependence of the central temperature on 
the mass. This can only be explained by a strong 
temperature dependence of the source of stellar 
energy, a fact which has not been sufficiently 
realized in the past. Theoretically, the central 
temperature increases somewhat with increasing 
mass of the star, more strongly with the mean 
molecular weight, and is practically independent 
of the chemical composition, i.e., of y and z. 

The radius of the star is larger for heavy stars 
and for high molecular weight. The density 
behaves, of course, in the opposite way. Both 
these results are in qualitative agreement with 
observation. The product of mass and density 
which occurs in Eq. (51) for the luminosity, is 
almost independent of the mass; therefore, for 
constant concentrations z, the luminosity is 
determined by the central temperature alone. 
Both radius and density are almost independent 
of the chemical composition, except insofar as it 
affects u. | 

The luminosity increases slightly faster than 
the fourth power of the mass and the sixth power 
of the mean molecular weight. This increase is 
considerably less than that usually given 
(M* 57-5) and agrees better with observation. 
The difference from the usual formula.is mainly 
due to the different dependence of the opacity on 
density and temperature; in fact, with the usual 
assumption (a=1, B=33), we obtain Mey, 
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_ The remaining difference is that the dependence 
on the radius is carried as a separate factor in the 
usual formula while we have expressed R in 
terms of M and yw. The observations show, for 
bright stars, an even slower increase than M?*; 
this seems to be due to the lower average molecu- 
lar weight (higher hydrogen content) of most 
very bright stars. It may be that these stars 
become unstable because of excessive energy 
production when their hydrogen content becomes 
too low (cf. §12).—The luminosity is inversely 
proportional to the concentration y of heavy 
elements (Russell mixture) because y determines 
the opacity. However, L is almost independent of 
the nitrogen concentration, as are R, T and p. 

All these considerations are valid for bright 
stars, down to about three magnitudes fainter 
than the sun. For fainter stars, with lower central 
temperatures, the proton combination H+H 
=D-+e+t should become more probable than the 
carbon-nitrogen reactions, because this reaction 
depends less on temperature. In discussing the 
energy production from H+H, we must take 
into account that (cf. §5) at low temperatures, 
this reaction leads only to He? rather than He‘, 
because of the slowness of the reaction He*+ He? 
= Be’. (Assumption B of §5 is made, viz. that He® 
is more stable than H’, and Li‘ unstable.) This 
causes a rather sudden drop in the energy 
evolution from H+H around 14-10° degrees 
(cf. Fig. 1), i.e., just below the temperature at 
which the proton combination becomes im- 
portant (~16-10° degrees, see Fig. 1). Therefore, 
the temperature exponent y stays fairly large 
(~13, cf. Fig. 2) down to about 13-10° degrees 
which corresponds to an energy production of 
about one percent of that. of the sun (five 
magnitudes fainter). For still fainter stars, y 
drops to very low values, reaching a minimum of 
-about 4.5 near 10’ degrees. 

The relations between central temperature, 
radius, luminosity and mass for this case (y = 4.5) 
are given in the last column of Table IX. The 
temperature is seen to depend much more 
strongly on the mass (as M®8) while the radius 
becomes almost independent of M and the 
density decreases with decreasing mass. The 
luminosity decreases faster with the mass (as M5) 
than for the bright stars. Unfortunately, little 
material is available for these fainter stars. This 
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is the more regrettable as it is rather important 
for nuclear physics to decide whether the H+H 
reaction is really as probable as assumed in this 
paper: There is some possibility that it is 
forbidden by selection rules (cf. reference 16, p. 
250) in which case it would be about 10° times 
less probable. Then the carbon-nitrogen reactions 
would furnish the energy even for faint stars, and 
the central temperature of these stars would not 
depend much on their mass. 

Figure 1 gives the energy production due to the 
proton combination (H+H) and to the carbon- 
nitrogen reactions (N+H) as a function of the 
central temperature, on a logarithmic scale. The 
great preponderance of H+H at low and N+-H 
at high temperatures is evident. The following 
assumptions were made: Hydrogen concentration 
35 percent, nitrogen 10 percent, central density 
p=100, average energy production = ¢ of central 
production for H+H, 5 for N++H reaction. He? 
is supposed to be more stable than H’, and Li‘ 
unstable (assumption B of §5). 


ENERGY 
PRODUCTION 





Fic, 2. The exponent y in the relation L~TY between 
luminosity L and temperature T, as a function of T. Solid 
curve: y for total energy production (logarithmic derivative 
of solid curve in Fig. 1). Dotted curves: y for stability 
against temperature changes (curve Ba for times less than 
14 months, Bb for more than 14 months). 
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Figure 2, solid curve, gives the exponent y in 
the total energy production. It is low at low 
temperatures (~ 4-5, hydrogen reaction) and has 
a minimum of 4.44 at 11 million degrees. Between 
11 and 14 million degrees, y rises steeply as the 
reaction Het-++He? = Be’ sets in (§5); then it falls 
again from 13 to 12 when this reaction reaches 
saturation. From T=15, the carbon-nitrogen 
reactions set in, causing a rise of y to a maximum 
of 17.5 at T=20, while at higher temperatures y 
decreases again as T~. 


§11. STABILITY AGAINST TEMPERATURE CHANGES 


Cowling has investigated the stability of 
stars against vibration. This stability is de- 
termined mostly by the ratio y of the specific 
heats at constant pressure and constant volume. 
If the radiation pressure in the star is negligible 
compared with the gas pressure (y=5/3) then 
the star will be stable for any value of the 
temperature exponent in the energy production 
(38), up to +450. Only for very heavy stars, for 
which the radiation pressure is comparable with 
the gas pressure, does the stability condition put 
a serious restriction on the temperature depend- 
ence of the energy production. According to our 
theory, the energy production is proportional to 
T!" (see below); according to Cowling, stability 
will then occur when y > 10/7. The corresponding 
ratio of radiation pressure to total pressure is 


Soy 
peg See 4 (60) 
3(7y—9) 

and the corresponding mass of the star, according 
to Eddington’s “standard model,” is 17/p? ©, 
where u is the average molecular weight. There- 
fore practically all the stars for which good 
observational data are available will be stable. 

It may be worth while to point out that the 
temperature exponent n for these stability con- 
siderations is not exactly equal to that for the 
energy production in equilibrium. A change of 
temperature gives rise to radial vibrations of the 
star whose period is of the order of days or, at 
most, a few years. On the other hand, when the 
temperature is raised or lowered, the equilibrium 


48 T, G. Cowling, Monthly Not. 94, 768 (1934); 96, 42 
(1935). 
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between carbon and nitrogen will be disturbed 
and it takes a time of the order of the lifetime of 
C! (~10° years) to restore equilibrium at the 
new temperature. Thus we must take the concen- 
trations of carbon and nitrogen corresponding to 
the original temperature. 

At T=2.-10’ degrees, the carbon reactions have 
a 7=50.6 (Table V) and therefore a temperature 
exponent yc = 16.2 (cf. 52); the nitrogen reactions 
have 7=56.3 and yn=18.1. The number of 
reactions per second is, in equilibrium, the same 
for each of the reactions in the chain (42). The 
energy evolution in the first three reactions 
together is 11.7 mMU, after subtracting 0.8; 
mMU for the neutrino emitted by N®. These 
three reactions are carbon reactions, the re- 
maining three, with an energy evolution of 
15.0 mMU, are nitrogen reactions. Thus the 
effective temperature exponent for stability 
problems becomes 


11.7y¥o+15.0yn 


y= =17.3 
26.7 


(60a) 


at 2-10" degrees. y is approximately proportional 
to T= (cf. 52, 6). 

For the proton combination, we have to 
distinguish the three possibilities discussed in 
§5. The simplest of these is assumption C. 


Assumption C: Lit unstable, H? more stable than He’ 


In this case, there are two “slow” processes, vtz. (a) the 
original reaction H+H and (b) the transformation of He’ 
into H? by electron capture (~2000 years). (a) depends on 
temperature approximately as T3, (b) is independent of T. 
The energy evolution up to the formation of He? is 7.2 
mMlU, from He? to Het 21.3 mMU. Thus 


_7.2-3.5+21.3-0_ 


0.9. 61C 
28.5 ere 


Y 


This would be a very slight dependence indeed, 


Assumption A: Li* stable 


According to Table V, the transformation of He? into 
Li‘ takes about one day. 

(a) For times shorter than one day, the reactions up to 
He?’ and from then on are independent of each other. The 
first group again has yH=3.5 and gives 7.2 mMU, the 
second group now gives only 8.3 mMU because the remain- 
ing 13.0 are lost to the neutrino from Li‘ (cf. §5) and has 
(Table V) yne=6.9. Therefore 


_7.2-3,54+8.3-6.9 


15.5 5.3. (61Aa) 


Y 
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(b) For times longer than one day, the proton combina- 
tion determines the whole chain of reactions so that 


y=3.5. (61Ab) 


Assumption B. Lit unstable, He? more stable 
than H? (most probable assumption) 


As was pointed out in §5, the reaction Het+ He’ 
= Be’ is so slow (3-10’ years, Table V) that the 
concentration of He? will remain unaffected by 
the temperature fluctuation. Be” has a mean life 
of 14 months at the center of the sun so that there 
are again two cases: 

(a) For times less than 14 months, there are 
three groups of reactions, (a) those up to He’, 
giving again 7.2 mMU with yu=3.5, (b) the 
reaction He++He*= Be’, giving (Table V) 1.6 
mMU with yoe=15.1, and (c) the electron 
capture by Be’, followed by Li’7+H =2Het-*. This 
last reaction does not depend on temperature 
(ype=0) and gives 18.6 mMU. Therefore the 
effective y becomes 


7.2:-3.5+1.6-15.1+18.6-0 
ye Se ene 


(61Ba) 
27.4 


(b) For times longer than 14 months, the reac- 
tion He*+ He’ = Be’ governs all the energy evolu- 
tion from He? on, so that 


7.2°3.5-4+20.2°15.1 


y= = 12.1 
27.4 


(61Bb) 


at 2-10’ degrees. 

At low temperatures, the reaction He*+ He? 
= Be’ stops altogether (cf. §5) so that then the y 
of the H+H reaction itself determines the radial 
stability. From 12 to 16-108 degrees we have a 
transition region in which the importance of the 
He‘+He? reaction (and the consequent ones) is 
reduced. l 

In Fig. 2, curves Ba and Bb, we have plotted 
the effective y for the radial stability, by taking 
into account both C+N and H+H reactions, 
and making the same assumptions about the 
concentrations of hydrogen and nitrogen as in 
Fig. 1 (cf. end of §10). Assumption B was made 
regarding the stability of Lit and He?; the curves 
Ba and Bb refer to the short time and long 
time formulas (61Ba) and (61Bb). At high 
temperatures, the two curves coincide because 
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then the proton combination is unimportant 
compared with the carbon-nitrogen reactions. 
The combined y is seen to reach a maximum of 17 
(for the long time curve). 


§12. STELLAR EVOLUTION“ 


We have shown that the concentrations of 
heavy nuclei (Russell mixture) and, therefore, 
also of nitrogen, cannot change appreciably 
during the life of a star. The only process that 
occurs is the transformation of hydrogen into 
helium, regardless of the detailed mechanism. 
The state of a star is thus described by the 
hydrogen concentration x, and by a fixed 
parameter, y, giving the concentration of Russell 
mixture. The rest, 1—x— y, is the helium concen- 
tration. Without loss of generality, we may fix 
the zero of time so that the helium concentration 
is zero. (Then the actual “birth”’ of the star may 
occur at t>0). 

It has been shown (Table IX) that the lumi- 
nosity depends on the chemical composition 
practically only*® through the mean molecular 
weight uw. This quantity is given by 


1/u=2x+3(1—x—y)+3y= (5/4) (x+a), (62) 
a=0.6—0.2y, (62a) 


taking for the molecular weights of hydrogen, 
helium and Russell mixture the values 4, 4/3 and 
2, respectively.! 

Now the rate of decrease of the hydrogen 
concentration is proportional to the luminosity, 
which we put proportional to u”. According to 
Table IX, n is about 6. Then | 


dx /dt ~—(x+a)-". (63) 
Integration gives 
(x+a)" t! =A (to— t), (64) 


where A is a constant depending on the mass and 
other characteristics of the star. Since x=1—y 
at t=0, we have 


Ato= (1.6—1.2y) "+1, (64a) 


44 Most of these considerations have already been given 
by G. Gamow, Phys. Rev. 54, 480(L) (1938). 

46 Except for the factor y7 which, however, does not 
change with time. 
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It is obvious from (63) and (64) that the 
hydrogen concentration decreases slowly at first, 
then more and more rapidly. E.g., when the 
concentration of heavy elements is y=4, the 
first half of the hydrogen in the star will be 
consumed in 87 percent of its life, the second half 
in the remaining 13 percent. If the concentration 
of Russell mixture is small, the result will be 
even more extreme: For y=0, it takes 92 percent 
of the life of the star to burn up the first half of 
the hydrogen. Consequently, very few stars will 
actually be found near the end of their lives even 
if the age of the stars is comparable with their 
total lifespan fo (cf. 64a). In reality, the lifespan 
of all stars, except the most brilliant ones, is long 
compared with the age of the universe as deduced 
from the red-shift (~2-10° years): E.g., for the 
sun, only one percent of the total mass trans- 
forms from hydrogen into helium every 10° years 
so that there would be only 2 percent He in the 
sun now, provided there was none “in the 
beginning.” The prospective future life of the sun 
should according to this be 12-10° years. 

It seems to us that this comparative youth of 
the stars is one important reason for the existence 
of a mass-luminosity relation—if the chemical 
composition, and especially the hydrogen con- 
tent, could vary absolutely at random we should 
find a greater variability of the luminosity for a 
given mass. 

It is very interesting to ask what will happen 
to a star when its hydrogen is almost exhausted. 
Then, obviously, the energy production can no 
longer keep pace with the requirements of equi- 
librium so that the star will begin to contract. 
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(This is, in fact, indicated by the factor g!/8 in 
Eq. (55) for the stellar radius; z is proportional to 
the hydrogen concentration.) Gravitational at- 
traction will then supply a large part of the 
energy. The contraction will continue until a new 
equilibrium is reached. For “light” stars of mass 
less than 6u7? sun masses (reference 1, p. 507), 
the electron gas in the star will become degenerate 
and a white dwarf will result. In the white dwarf 
state, the necessary energy production is ex- 
tremely small so that such a star will have an 
almost unlimited life. This evolution was already 
suggested by Strémgren.! 

For heavy stars, it seems that the contraction 
can only stop when a neutron core is formed. The 
difficulties encountered with such a core*® may 
not be insuperable in our case because most of 
the hydrogen has already been transformed into 
heavier and more stable elements so that the 
energy evolution at the surface of the core will 
be by gravitation rather than by nuclear reactions. 
However, these questions obviously require much 
further investigation. 

These investigations originated at the Fourth 
Washington Conference on Theoretical Physics, 
held in March, 1938 by the George Washington 
University and the Department of Terrestrial 
Magnetism. The author is indebted to Professors 
Strömgren and Chandrasekhar for information 
on the astrophysical data and literature, to 
Professors Teller and Gamow for discussions, and 
to Professor Konopinski for a critical revision of 
the manuscript. 


46 G., Gamow and E. Teller, Phys. Rev. 53, 929(A), 


608(L) (1938). 


